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Abstract 



This article extends the previous paper in "M.W. Yuen, Stabilities for Euler-Poisson Equa- 
tions in Some Special Dimensions, J. Math. Anal. Appl. 344 (2008), no. 1, 145-156.", from 
the Euler-Poisson equations for attractive forces to the repulsive ones in R N (N > 2). The 
similar stabilities of the system are studied. Additionally, we explain that it is impossible to 
have the density collapsing solutions with compact support to the system with repulsive forces 
for 7 > 1. 
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1 Introduction 

O 

The semi-conductor models can be formulated by the isentropic Euler-Poisson equation with re- 
pulsive forces in the following form: 



p t +V ■ {pu) =0, 

(ptt) t +V • (pu <g> u) + VP + Ppu =pV$, (1) 



I A$(t,x)=a(N)p, 
where a(N) is a constant related to the unit ball in R N : a(l) = 2; a(2) = 2ir; For N > 3, 

n N/2 

a(N) = iV(jV - 2)7(AQ = JV(JV - 2) - , (2) 

where V(iV) is the volume of the unit ball in R N and T is the Gamma function. As usual, 
p = p(t,x) and u = u(t,x) G are the density and the velocity respectively. P = P(p) is the 
pressure and /3 > is the frictional damping constant. In the above system, the self-repulsive 
potential field $ = $(t, x) is determined by the density p through the Poisson equation. The 
equations <jTJ) i and (fljlo are the compressible Euler equations with forcing term. The equation (JT|) 3 
is the Poisson equation through which the potential with repulsive forces is determined by the 
density distribution of the electrons themselves. Thus, we called the system ([T]) the Euler-Poisson 
equations with repulsive forces. In this case, the equations can be viewed as a semiconductor 
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model. See [2], [5] for a detail about the system. For some fixed K > 0, we have a 7-law on the 
pressure P(p), i.e. 

P{p) = Kp<=^, (3) 

7 

for K > 0, which is a common hypothesis. When K — 0, the pressureless system can be used as 
models in plasma physics [T]. The constant 7 = Cp/c v > 1, where cp and c„ are the specific heats 
per unit mass under constant pressure and constant volume respectively, is the ratio of the specific 
heats. Additionally, the fluid is called isothermal if 7 = 1. 
The Poisson equation can be solved as 

§(t,x)=[ G(x-y)p(t,y)dy, (4) 

JR N 

where G is the Green's function for the Poisson equation in the ./V-dimensional spaces defined by 

M, N = l; 

G(x) = { log |a;|, N = 2; (5) 

1 /2 

If we seek solutions in radial symmetry with the radius r = f^Z^Li xf\ , the Poisson equation 
(jT])3 is transformed to 

r"" 1 ^ (t, x) + (N- 1) r N - 2 $ r =a (N) pr N ~\ (6) 

^=^5 f p(t,s)s N -'ds. (7) 
r Jo 

We can seek the radial symmetry solutions 

pit, x) = pit, r) and u = -V(t, r) =: -V. (8) 

r r 

By standard computation, the Euler-Poisson equations in radial symmetry can be written in the 
following form: 

N — 1 

Pt + V Pr + pV r + —^—PV = 0, ^ 

p(V t + VV r ) + P r (p) = p$ r (p). 

Perthame discovered the blowup results for 3-dimensional pressureless system with repulsive forces 
[S]. In short, all the results above rely on the solutions with radial symmetry: 

V t + ^ r =^ /; p{t, s)s N -Hs. (10) 

And the Emden ordinary differential equations were deduced on the boundary point of the solutions 
with compact support: 

= i? ^T' R (°> R o) = Ro>0, R(Q, Ro) = 0, (11) 

where := V and M is the mass of the solutions, along the characteristic curve. They showed 
the blowup results for the C 1 solutions of the system ©. 

Very recently, Yuen showed in [T3] that the classical non-trivial solutions (p, V) for the Euler or 
Euler-Poisson equations with repulsive forces, in radial symmetry, with compact support in [0, R], 
where R is a positive constant (which is V(t,0) = 0;p(t,r) = 0, V(t,r) = for r > R) and the 
initial velocity such that: 

H = [ V dr>0, (12) 
Jo 
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blow up on or before the finite time T — 2R/Hq, for pressureless fluids (K = 0) or 7 > 1. 

The system with attractive forces were studied in [3], [B], [7], jS], |12).and [13]. This article extends 

the previous paper in [13] , from the Euler-Poisson equations for attractive forces with or without 

frictional damping to the repulsive ones in R N (N > 2). The similar stabilities of the system are 

studied. 

In the last section, we exclude the possibility of collapsing solutions for this system. The 
non-existence of collapsing solutions can be shown by the simple argument for the energy function: 

E (t) = In (hP M 2 + ^r P " dx > for 7 > 1- (13) 

Theorem 1 For the classical solutions with compact support of the Euler-Poisson equations with 
repulsive force, (Qp, in R N (N > 2) with 7 > 1 or without pressure (K — 0), there is no collapsing 
phenomenon where part of the density pit, x) collapses to a point. 



2 Stabilities 



In this section, we study the stabilities of the Euler-Poisson equations with repulsive forces, ([T]), 
in R N (N > 2). The total energy can be defined by, 



E(t) = In(y\ u \ 2 



VM 2 + ip- ipfcW for 7 >l, 



5/0$) dx, for without pressure. 



(14) 



For the system, we have the energy estimate: 



Lemma 2 For the Euler-Poisson equations, (J\), suppose the solutions (p, u) have compact support 
in CI. We have, 

|2 

\u 

in 



E(t) = -P / p\u\ dx < 0. 
Jn 

Initially, Sideris used the second inertia function 



(15) 



H(t) 



p(t, x) \x\ dx, 



(16) 



to study instability results for the Euler equations [TU]. After that in [5J, the instability result of 
the Euler-Poisson equations with attractive forces, in radial symmetry, was obtained for 7 > 4/3 
and N = 3. For the corresponding cases in R N , with non-radial symmetry, were studied [3J, [13] . 
By the standard computation for energy method, it is clear to have the following lemma: 

Lemma 3 Consider {p, u) is a solution with compact support in CI for the Euler-Poisson equations, 
(QJj with /3 = 0. We have 



H(t) = 2 L (p\u\ 2 +NP)dx 



N-2 



dx, for N > 3; 



H(t) = 2 J n (p \u\ 2 + 2P)dx + M 2 , for N = 2. 
By applying the above lemma, we could have: 



(17) 



Theorem 4 Suppose (p, u) is a global classical solution in the Euler-Poisson equations, (0J) with 
7 > 1 , without frictional damping (/5 = 0) . We have 
(l)for N > 3, 



liminfM > 

t->oo t 



inf(2,iV(7- \),N -2)E 
M 



-1 1/2 
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(2)for N = 2, 



(S)for N>2, 



lim inf ^ >J^; 
t^oo f - V 2M' 



lim inf^> 

t— >oc f 



NKM"!- 1 



(7-1) 



1/2 



wi£/i R(t) = max l£ f!( t ) {M}- Here 



M = / p(t, x)dx, 
Jn 



(19) 



(20) 



(21) 



is the total mass which is constant for any classical solution and is the fixed volume oftt. 
Proof. (l)For N > 3, we have the positive energy function E > 0. We can get from Lemma [31 

N -2 





p\u\ 2 + NP 







dx 



That is 



H(t) > H(0)+H(0)t + inf (2, N(j — 1),N — 2)Et 
On the other hand, we obtain, 



p$dx \ > 2inf(2,7V( 7 - 1),N - 2)E. (22) 

(23) 



That is, 



H(0)+H(0)t + inf(2, N(j — 1),N — 2)Et 2 < H (t) < R{t) 2 M. 



0( 1 -)+mi(2,N( 1 -l),N-2)E< R( ^, 



lim inf M> 



inf(2,iV(7-l),iV-2)J5 



1/2 



For N = 2, we have 



#(*) = 2 / (p |< + 2P)dx + Af 2 > M , 
~ i 2 + C a t + Ci < (t) < i?(t) 2 M, 



For AT > 2, we obtain 



and 



hminfM>J^. 
t^oo f - V 2M 



l/ 7 



,(7-l)/7 





p\u\ 2 + NP 







d.r - — 2 I p<S>dx } > 2 I NPdx. 

Q 



From the inequality ([30]) . it is clear to have 

Hit^NKlfl] 1 ^ AP> 0, 



H(0)+H{0)t + NK MV< H(t) < R(t) M, 
t t z 



(24) 

(25) 
(26) 

(27) 
(28) 
(29) 

(30) 
(31) 

(32) 

(33) 
(34) 
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This gives 

r R(t) \NKM~t- 11 ' " 
lim inf > -, — T7— 



(35) 



The proof is completed. ■ 

By applying the similar method, in the above section to the system, ([T]), with frictional damping 
constant (/3 > 0), the below theorem is followed: 

Theorem 5 Suppose (p, u) is a global classical solution with compact support in the system, (QJ), 
with frictional damping (/3 > 0). We have for N > 2, 

R(t) ( 2BNKM~<- 1 \ 1/2 , 
lim inf — - ( , 36 

with R(t) = max xe n(t) {M}- 
Proof. For > 2, we get, 

P(i) = / 2xpudx, (37) 

and 

#"(*) = 2 [ x [-V • (pu (g) u) - VP + - ppu] dx. (38) 
Additionally, it can be arranged as the following: 

H(t) = 2 I [x - V • (pu ® u) - VP + pV$] dx - 0H(t), (39) 

H(t) + ^H(t) = 2 ij p \u\ 2 + NP dx - J p$rfa;| > 2 J NPdx > 2NK | O | 1 ^ M 7 > 0. 

(40) 

Therefore, we are able to obtain the inequality, 

C 3 +Cie~ pt +20NK |r2| 1_7 NPt< H(t) < R(t) 2 M, (41) 

O(-) + l^l 1 ' 7 M 7 < gWJ^ , (42) 

This gives 

,. . ,P(t) ( 2(3NKM~<- 1 \ 1/2 

Km inf -p^> — — , . 43 

The proof is completed. ■ 



3 Non-existence of Collapsing Solution 

In this section, we explain the idea that there is no possibility to have a density collapsing solution, 
with compact support for the Euler-Poisson equations with repulsive forces: 
We restate their energy for 7 > 1 or the pressureless fluids is: 

E{t) = f (\p M 2 + ~^lP - \p®) dx>- f l -p<S>dx. (44) 

When a <5-shock exists for the density function p(t,x), the potential energy function, with N >3, 
becomes to be infinite: 

- f p^dx= f p{t,x) f p( " t,V }_ 2 dydx = lim f S(t,x) f ^j^-dydx = 00. (45) 
Jn Jn Jn \x — y\ c->o+ J n Jq e 
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With N = 2, the situation is similar: 

— / pQdx = — / S(t,x) / S(t, y) In \x — y\ dydx = oo. (46) 
Jn Jn Jn 

Therefore, the total energy of the <5-shock density solutions must be infinite. However, for the 
classical solutions with compact support, the initial energy is finite. By the energy estimate in 
Lemma [2 we have, 

E(t) < E(0). (47) 

If the total energy is finite, it is impossible to obtain the density collapsing solutions. It is clear to 
have Theorem [TJ 
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